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We study preheating in theories where the infiaton couples derivatively to scalar and gauge fields. 
Such couplings may dominate in natural models of inflation, in which the flatness of the infiaton 
potential is related to an approximate shift symmetry of the infiaton. We compare our results 
with previously studied models with non-derivative couplings. For sufficiently heavy scalar matter, 
parametric resonance is ineffective in reheating the universe, because the couplings of the infiaton 
to matter are very weak. If scalar matter fields are light, derivative couplings lead to a mild long- 
wavelength instability that drives matter fields to non-zero expectation values. In this case however, 
long-wavelength fluctuations of the light scalar are produced during inflation, leading to a host of 
cosmological problems. In contrast, axion-like couplings of the infiaton to a gauge field do not 
lead to production of long-wavelength fluctuations during inflation. However, again because of the 
weakness of the couplings to the infiaton, parametric resonance is not effective in producing gauge 
field quanta. 



I. INTRODUCTION 

During a sufficiently long epoch of inflation, matter 
is diluted away by the quasi-exponential expansion of 
the universe. The re-population of the universe with 
radiation after the end of inflation is then typically 
achieved by the decay of the inflaton into ordinary mat- 
ter particles. The perturbative theory of this process — 
reheating |TJ[2] — was refined by the discovery [31I11[S] that 
in many models the dynamics would proceed through a 
stage of parametric resonance, leading to the extremely 
efficient decay of the inflaton into a far from equilibrium 
distribution of matter. The understanding of this prelim- 
inary stage of reheating — preheating — has since been de- 
veloped by many authors [ilTllHlElIinillllliaiSlIIlIIS]- 

Most work on preheating has focused on models with 
direct, non-derivative, couplings to matter. Little is 
known about the strength of these couplings, but the- 
oretical arguments suggest that they have to be rather 
weak. Indeed, in order for the inflaton to drive a phe- 
nomenologically acceptable stage of inflation, its poten- 
tial has to be extremely "flat". A variety of different 
ways to stabilize these flat potentials have been stud- 
ied [IH], but perhaps the most compelling idea that has 
emerged to date is that inflation is driven by a pseudo- 
Nambu-Goldstone boson [VT^ ■ In the simplest realization 
of this idea, the inflaton sector is invariant under an ap- 
proximate global U{1) symmetry that shifts the inflaton 
field by a constant. Because only a constant potential is 
invariant under shifts of the infiaton, this approximate 
symmetry guarantees that deviations from fiatness are 
small. But at the same time, because direct couplings 
between the infiaton and matter violate the shift sym- 
metry, they are are also expected to be extremely weak, 
and perhaps even negligible. 
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On the other hand, derivative couplings of the infiaton 
to matter do satisfy the shift symmetry, and so there is 
no reason for them to be particularly weak. It is there- 
fore entirely possible that derivative couplings could be 
more important than non-derivative ones during reheat- 
ing. Studies of preheating in models with derivative cou- 
plings 1^ HI] have been less extensive than those of di- 
rect couplings, and therefore our aim here is to carry 
out one such study: examining preheating in these mod- 
els and determining whether they lead to a qualitatively 
different picture of the end of inflation and the onset of 
reheating. 

In the next section we will review the standard ideas 
of preheating in models with direct couplings between 
the inflaton and matter. In section UTTl we then motivate 
the study of derivatively coupled models and carry out 
the related preheating calculations, including couplings 
to both scalar and gauge flelds, before concluding. 



II. PARAMETRIC RESONANCE IN 
CANONICAL MODELS 

We begin by reviewing the basic results about preheat- 
ing in a simple model with a direct coupling between the 
inflaton and matter flelds. For detailed methods and re- 
sults regarding the physics of reheating and preheating, 
we refer the reader to [15,, 20J . 

For simplicity we assume inflation is driven by a single 
real scalar field (f>, slow-rolling down a quadratic effec- 
tive potential, and coupled to a massive real scalar field 
X (representing matter fields) through a quartic interac- 
tion, 

(1) 

Throughout we treat the infiaton field as spa- 
tially homogeneous and work in a spatially fiat 
Friedmann-Robertson- Walker universe with metric 
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We expand the field x i^i terms of annihilation and 
creation operators 



(Pk 
(27r)3/^ 



(afcXfeWe**'-" + h.c.). (2) 



Then the mode functions Xfc and the homogeneous infla- 
ton (/> satisfy the field equations 



Xk + 3ffxfe 



2 J,2 



Xfc=0, (3) 
(4) 



where the Hubble parameter H(t) = a/a is determined 
by the Friedmann equation 



1 



2 j,2 



(5) 



and where we have denoted a time derivative by an over- 
dot. Note that we have neglected the effects of the x field 
and its interactions with on the expansion rate of the 
universe, since we are interested in an epoch at which the 
inflaton is dominant, so that <C p^. Similarly, we will 
henceforth also neglect the back-reaction term g^x^cj) ™ 
equation Q. 

Slow-roll inflation occurs when <g; m^(j) f^'^d 
(j) <C 'iHcf). During slow-roll, the interaction between 4> 
and X in equation (|3| simply increases the effective mass 
of the field x. 



2 2 , 2i2 



(6) 



There are several arguments that suggest that this mass 
should not be much lower than the Hubble scale. First, 
if TTtcff <^ H, a, non-zero value of x remains frozen during 
inflation, so there is no reason to assume that x lies at 
the minimum of its potential. At the same time, if the 
effective mass of x is smaller than the Hubble scale, long- 
wavelength fluctuations of x created during the inflation- 
ary epoch generically lead to isocurvature perturbations. 
And finally, if rrics ^ H, the same long- wavelength fluc- 
tuations will drive x away from zero, even if x = ini- 
tially, which could eventually lead to a second stage of 
inflation driven by the fleld x [21^- If "m^ is negligible, 
the requirement H < rrios, in combination with equation 
([5]), leads to the bound 



> ^ 
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(7) 



The same bound also guarantees that the effective mass 
of Xi equation ^ , is adiabatically constant during infla- 
tion, even if is negligible. Hence, the coupling be- 
tween the inflaton and any other light scalar cannot be 
too weak. 

Once inflation ends, the inflaton field begins oscillat- 
ing around the minimum of its effective potential. Dur- 
ing this oscillation stage, the expansion of the universe 



mimics a matter domination phase. In particular, soon 
after the end of inflation, the fleld evolution follows the 
equation 



(bit) « 0e 



(8) 



where 



Mpi denotes the value of the inflaton at the 



beginning of the oscillating stage^, at i « w^^^. Because 
of the rapid oscillations, the effective mass of x ceases 
being adiabatically constant for sufficiently light fields. 
One may then ask whether this non-adiabatic evolution 
leads to particle production. The answer is affirmative 
for appropriate couplings between matter and the infia- 
ton. More importantly, parametric resonance can excite 
X very efficiently. 

If we neglect the expansion of the universe and the 
back-reaction from x, equations Q and ^ can be com- 
bined into the Mathieu equation. 



Xfc + [<5-2ecos(2t)]xfc = 0, 



(9) 



where we have defined a dimensionless time t — s- t, 
and used 



6 = 



k 

m0 



- 2e and e — 



2m^ 



(10) 



Solutions of the Mathieu equation are known to exhibit 
parametric resonance — resonance for certain values of the 
dimensionless parameters S and e. In the S-e plane, these 
resonant solutions form band-like patterns called insta- 
bility bands (see ^2]). Along these unstable solutions, 
the mode functions grow exponentially. 



Xk oc exp{^j,m^t), 



(11) 



where the characteristic exponent /i depends on 6 and e. 

Although the Mathieu equation cannot be solved ex- 
actly for all values of d and e, using perturbative tech- 
niques in the regime of small e, one can show that the 
tips of the instability bands occur at 6 — where n is an 
integer [22l [S^ . A numerical plot of the bands is shown 
in Fig. [T] Extrapolating the tips of the plotted instabil- 
ity bands down to the e = axis shows agreement with 
the above perturbative results. In the theory we consider 
here however, equation ([t]) implies e ^ 1. Therefore, pre- 
heating occurs in the broad-resonanace regime, which, as 
we shall see, is quite different from its narrow-resonance 
counterpart at e < 1. In Fig. [2] we give examples of the 
behavior of unstable and stable modes. Similar results 
and a more in-depth discussion can be found in |15j . 

Since we have neglected the expansion of the uni- 
verse, the parameters of the Mathieu equation are time- 
independent. A mode that lies within one of the insta- 
bility bands therefore continues to grow indefinitely. Be- 
cause the effective particle number density Uk is related 



^ More precisely, the value of the inflaton at the beginning of the 
oscillating phase is actually smaller, <f'c ~ O.SMpi. 
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FIG. 1: Instability diagram for solutions of the Mathieu equa- 
tion in the S-e plane. In this diagram, shaded regions indi- 
cate resonant, or unstable, solutions of the Mathieu equation. 
Unshaded regions indicate stable solutions. We also show 



the evolution of the parameters in equation ( 14 1 for differ- 
ent comoving modes and fields: k — 0,m^ = (continuous), 
k = Q, 7^ (dashed) and k ^ 0, 7^ (dot-dashed). 



to the mode functions by 



nk 



1 

2u;i 



(IXfc 



2' 



(12) 



where = m^g -t- (k/a)"^, parametric resonance leads to 
unbounded growth of particle number, as shown in Fig. [3] 
Of course this unending particle production is simply a 
consequence of neglecting any dissipative effects on the 
inflaton (the back-reaction of matter and the expansion 
of the universe) . 

When the expansion of the universe is included, the 
inflaton undergoes damped oscillations resulting in a loss 
of energy with which to drive parametric resonance in 
the matter flelds. But on the other hand, the expansion 
of the universe redshifts all momenta, which causes the 
energy of a mode to decrease as the universe expands. 
It is possible to see both of these effects analytically by 
making the approximation a oc t^^^ and substituting into 
equations ^ and (|3|. Then, the rescaled field 



Xk = Xka 



3/2 



(13) 



satisfies the Mathieu equation ([9]), but now the parame- 
ters S and e become time-dependent 



S = 



m0 



2e and e = 



\2m^t ) 



(14) 

Even though b and e are time-dependent, their values 
are nearly constant during the oscillating phase. For in- 
stance, the relative change in e during an oscillation of 



X, 




FIG. 2: Top panel: Time evolution of an unstable mode when 
expansion of the universe is neglected. Here = 0, fc = m^, 
and fl^ = 1 X 10"^ 

Bottom panel: Time evolution of a stable mode when expan- 
sion of the universe is neglected. Here = 0, fc = Sm^, 
and = 1 X 10~^. Here, and in all remaining figures, time is 
measured in units of 27r/m^ which corresponds to the number 
of oscillations in S from the end of inflation. 
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(15) 



Hence, at late times, t k, H ^ ^ we can think of 
the parameters 8 and e as being (locally) constant. The 
same assumption of late times has been actually made in 
the derivation of ( 14 1 . 

The time-dependence in b and e will cause comoving 
modes to migrate in the b-t plane. Thus, when the ex- 
pansion of the universe is taken into account, particle 
production in a given /c-mode may take place for a brief 
time interval as this mode passes in and then out of an 
instability band. We also plot the trajectories of different 
modes in the b-e plane in Fig. [T] The longer a partic- 
ular mode remains inside an instability band, the more 
efficient is the production of the corresponding particles. 
It is clear from the figure that parametric resonance is 
always more effective for longer wavelengths and lighter 



FIG. 3: Comoving number density of x-particles in an un- 
stable mode when expansion of the universe is neglected. 
Here, as in the top panel of Fig. [2] — 0, k — m^, and 
= 1 X 10"^ 
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fields [Hj. If the initial value of e is small, the trajec- 
tory of, say, the fc = mode only crosses the narrow tip 
of the first instability band; this is the narrow-resonance 
regime, which is not particularly effective in producing 
particles. On the other hand, if the initial value of e is 
sufficiently high, comoving modes remain inside several 
instability bands a significant fraction of the time. We il- 
lustrate the behavior of unstable modes in an expanding 
background in Fig. [4j Note that parametric resonance 
can excite modes of x with energies much higher than 
the infiaton mass if e is sufficiently big. 

Other models with canonical couplings between the in- 
ffaton and matter can be similarly examined. For ex- 
ample, one can include trilinear terms proportional to 
that would arise naturally in theories in which (j) 
has a non- vanishing vacuum expectation value. Although 
the quantitative results differ, similar qualitative results 
emerge. What is less explored is whether models with 
non-canonical couplings between the inffaton and matter 
yield similar results. In the next section we examine such 
couplings. 
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FIG. 4: Top panel: Time evolution of an unstable mode 
when expansion of the universe is included. Here m-^ = 0, 
k — 10~^ X m^, and = 1 x 10~^. The resonance is less 
explosive than when expansion is neglected, and it only lasts 
briefly as the mode passes through an instability band. After 
the mode passes through the band, it begins to decay due to 
the expanding background. 

Bottom panel: Comoving number density of x-particles in the 
same mode as top panel. Here it is clearly seen that the mode 
passes through two instability bands (the first between t ~ 8 
and t ~ 12, the second between t ~ 15 and t ~ 30). While 
it passes through each band the comoving number density 
increases exponentially. 



III. REHEATING IN MODELS WITH 
DERIVATIVE COUPLINGS 

A. Derivative Couplings in Inflationary Scenarios 

Phenomenologically viable inffationary models require 
very "ffat" potentials. A simple example is pro- 
vided by our chaotic model with quadratic potential 
V{cf)) — |m^0^, in which sufficient ffatness implies that 
the mass of the field must be m « 10~^Mpi. If the inffa- 
ton did not couple to matter, this value would not pose 
any particular problem, as no quantum corrections (other 
than gravitational) could drive m to higher values. How- 



ever, the inffaton must couple to matter for reheating to 
occur, and these couplings will generically be responsi- 
ble for large quantum corrections. In the absence of any 
symmetry, a light inffaton hence requires ffne tuning. 

To be more speciffc consider an inffationary model with 
bare Lagrangian given by equation ([T]). To lowest order 
in the effective mass of receives a correction given 
by the Feynman diagram in Fig. [5] The loop integral is 
quadratically divergent, so the corrected mass squared is 



9^ 
167r2' 



(16) 



where A is the ultraviolet cut-off of the theory, which is 
expected to be Planckian. There are hence two different 
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ways to obtain m « 10~^Mpi: (i) The coupling of the 
inflaton to matter is strong, g « 1, and there has to be a 
cancellation between the bare mass and the quantum cor- 
rections of 1 part in 10^, or (ii) the inflaton is very light 
and it couples very weakly to matter. Case (i) requires 
fine tuning; as we shall see, case (ii) is realized in models 
where the inflaton is a Pseudo-Goldstone boson. Inci- 
dentally, let us point out that if quantum corrections to 
the mass of the inflaton are to be small, g^A^ / 
the bound in equation Q has to be violated. 



iT:^ < ml, 




FIG. 5: One-loop quantum correction to the mass of the in- 
flaton. A quartic coupling between the inflaton (p and matter 
X leads to a quadratically divergent mass. 

Because of the fine tuning required in case (i), sig- 
nificant efforts have been devoted to develop inflation- 
ary models in which the inflaton is naturally light |16| . 
As in many standard model extensions, these models re- 
quire additional symmetries. Perhaps the most promising 
candidates are models in which the inflaton is a pseudo- 
Nambu-Goldstone boson, the Goldstone boson of an ap- 
proximate, spontaneously broken symmetry |17j . In the 
limit that the symmetry is exact, the Goldstone boson is 
massless, and the inflaton potential is exactly flat. Small 
symmetry breaking terms then give a correspondingly 
small mass to the inflaton, and the remaining approxi- 
mate symmetry guarantees that quantum corrections re- 
main small. ^ 

As we have discussed above, if the inflaton is to de- 
cay into matter flelds, it is, of course, important that it 
couple to matter. In the context of natural models, this 
generically implies that matter fields must be charged 
under the spontaneously broken symmetry. To begin our 
exploration of reheating in natural models of infiation, 
let us therefore consider the simplest model that satisfies 
both requirements: a theory with two charged complex 
fields <i> and invariant under a spontaneously broken 
(approximate) U{1) symmetry, 

(F2 



$*$)2 



(17) 



This Lagrangian is invariant under the U{1) transforma- 



tion $ 



e *'^x, where is a complex number. 



^ It is sometimes argued that quantum gravity effects, such as the 
creation of virtual black holes, explicitly break this symmetry 
|25l . We shall ignore this possibility here. 



This phase will later play the role of the inflaton, while 
X will play the role of matter. For simplicity, we also 
assume that the theory has an additional Z2 symmetry 
X -X- 

The U{1) symmetry is spontaneously broken, because 
in the vacuum the field $ has a non- vanishing expectation 
value, (|<f>|) = F. As we shall argue below, in order 
for the infiaton potential to be sufficiently fiat, we need 
F w Afpi. To study the low-energy excitations around 
the vacuum, we redefine 



^ ^ F exp 



X Xexp 



F 



(18) 



and ignore "radial" excitations around the minimum 
of the potential because they are extremely heavy — 
with masses of order F. Substituting these expressions 



into ( 17 1 we find 



'm\x*X 



F2 



(a^0)(9^0)lx' 



- pid^(^){x*d^x-xd''x*) 



(19) 



Note that this automatically generates derivative cou- 
plings between the inflaton </> and the fleld x (which could 
be redeflned away only if A = 0). These couplings are 
model-independent, in the sense that they originate from 
the kinetic term of the complex field x ^^nd the field re- 
definition ( 18 ). 

For to be a viable infiaton, it is necessary to generate 
a potential for it. To this end, we introduce terms in 
the Lagrangian that explicitly break the J7(l) symmetry. 
Because the potential has to vanish in the limit of exact 
symmetry, these terms will generically lead to an inflaton 
potential of the form 



F(0) « 



1 - cos I — 



(20) 



where /i is a parameter with dimensions of mass that 
characterizes the strength of the symmetry breaking 
terms. Note that the potential has to be periodic, since 
(j)/F is a phase. For (j)/F < 1, we can assume that the 
potential is quadratic, with squared mass = fi'^/F'^, 
which is the form we shall consider in the following. 

It is convenient to work with real fields, rather than 
complex ones. Defining 



X = Xr + iXi 



(21) 



where Xr and Xi Etre, respectively, the real and imaginary 
parts of X, the Lagrangian ( 19 1 becomes 



C 



(9^0)9^0 - {d,Xr)d''Xr - {^^X^)^^Xi 



{ml 
2 



2XF^)xl - {ml - 2\F^)xl 



(22) 



+ pid^'P) {Xr^>^X^ ~ X^^>^Xr) " ^ (5^</') (5^0) (x. + X 
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Therefore, derivative couplings between the inflaton and 
matter are quite generic and come in two classes: cubic 
and quartic. Note that in order to avoid a tachyonic 
instability, we must assume ml > 2XF. It is also easy 



to verify that the potential (20 1 will satisfy the slow- 
roll conditions only for F > Mp\. It is unclear whether 
we can trust an effective field theory description with 
F ^ Mpi, so we shall implicitly assume F « Alp. Since 
Mpi denotes the reduced Planck mass, we could push this 
choice a bit higher. 



B. Parametric Resonance 

In section |ll] we described how the inflaton may de- 
cay by resonantly exciting modes of its decay products. 
This process is non-perturbative in nature and requires 
us to solve the equations of motion of the matter fields x 
during the oscillating stage of the inflaton at the end of 
inflation. We now wish to examine the extent to which 
this phenomenon may occur in models with derivative 
coupHngs, such as (19). 



Because, by assumption, the fields x have a vanishing 
expectation value at the end of inflation, the cubic terms 
in the interaction will have no effect on their equation 
of motion. Thus, we may concentrate on the quartic 
interaction. 



C 



1 



1 



2 2 



1 

F^ 



(23) 



The resulting equations of motion are 



1 + 2 



F 



^{d^cl,)d^<l> + mlx, (24) 



with the mode equation for x therefore given by 

"2 „ " 



Xk + iHxk 



F2 



Xk = 0. (25) 



The form of this equation already hints that inflation 
and preheating in derivatively coupled models may be 
qualitatively different from the canonical case considered 
in section |TT] This is because the field x has an effective 
squared mass given by 



'off 



(26) 



Although (f)^ /F"^ is small during slow-roll inflation, the 
correction term becomes large when the inflaton leaves 
the slow-roll regime. Therefore, for sufficiently small val- 
ues of my, the effective mass of the field x might become 




FIG. 6: Instability diagram for solutions of the Mathieu equa- 
tion in the 5-e plane. We superimpose the evolution of the 



parameters in equation ( 28 1 for different comoving modes and 
fields: k = G,m^ = (continuous), k = Gyfriy^ 7^ (dashed) 
and k 7^ 0,771,^ 7^ (dot-dashed). 



tachyonic. Clearly, such an instability is absent as long 
as the effective mass is positive. 



> 



V2m0Mpi 



F 



(27) 



Interestingly, this bound is very similar to the require- 
ment that the mass of x be larger than the Hubble scale 
during inflation. As a matter of fact, since F « Mpi, 
the bound (27 1 simply reads > m^. During slow-roll 



H K, 0/Afpi, and because (j) is of order Mpi during 
the period of inflation accessible to present observations, 
equation (27 1 also follows from requiring > H. It 



turns out that the tachyonic instability associated with 
the violation of the bound (27 1 is actually harmless, and 



it is just the production of long-wavelength fluctuations 
of the field x what renders light scalar fields undesirable. 
To further study parametric resonance in this class of 



models, it will be instructive to recast equation (25 1 into 



the form of the Mathieu equation. Including the expan- 
sion of the universe but neglecting back-reaction from x, 
we find that the rescaled variable ( 13 1 satisfies the Math- 
ieu equation ([9|, with 



/ k 



2e and e 



Ft 



m 

The diagram in Fig. [6] shows different trajectories of 
a mode in the 6-e plane. It is again apparent from the 
figure that for a given initial value of e preheating is in 
general more efficient for longer wavelength modes and 
lighter fields. A crucial property of derivatively coupled 
models is that e is at most of order one, because (j)^ ~ 
Mpi w F. Hence, any resonance proceeds close to the 
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end of the instability bands, and is clearly ineffective if 
fields are heavy, ((5 » 1). On the other hand. 



if condition (27) is violated, long-wavelength modes will 



remain inside the instability band as long as the effective 
mass is negative, which could lead to a significant growth 
of the homogeneous component of %. 

In order to estimate the efficiency of preheating in the 
tachyonic regime, let us consider the evolution of the 
k — Q mode of a light field, <^ m^. As we argued 
above, this is the case where parametric resonance is 
strongest. After the end of inflation, the mode fc = 
finds itself in the instability band at 5 < 0. To lowest or- 
der in e, the characteristic exponent is simply n = y/—8. 
Hence, for adiabatic changes in /i, we expect the solu- 
tions to grow as exp(J dt). However, this is only a 
good approximation if 0c/^ 2> 1. Indeed the adiabatic- 
ity parameter is 



F 



(29) 



which is of order one instead. Therefore, we shall solve 
equation (|9| directly by approximating 



5^- 



Ft 



e = 0. 



(30) 



The solution is a linear combination of powers in t. Keep- 
ing the growing mode, and returning to the original vari- 
able, we find 



where 



P- 



p2 ' 



(31) 




FIG. 7: The zero mode of x for "T'x = lO^^m^. During the 
initial stages, the effective mass is tachyonic, and the mode 
grows as it remains in the stability band at 5 < 0. Once the 
mode exits the tachyonic regime, the mode function decays. 
Note that long-wavelength excitations do not admit a particie 
interpretation. 



with vector fields as the only remaining alternative. Inci- 
dentally, vector fields arc particularly appealing as decay 
products of the inflaton for two reasons: Their masses 
are protected by gauge symmetries, so they are natu- 
rally light, and the conformal nature of their couplings to 
gravity prevents the production of long-wavelength fluc- 
tuations of these fields during inflation. 

Consider therefore an effective theory in which a light 
pseudo-scalar (j) (the inflaton) couples derivatively to a 
U{1) gauge field (matter). 



and Xc denotes the value of x the end of inflation. 
Clearly, as seen in Fig. ([t]), for (fyc/F w 1 this ampli- 
fication is very modest. Hence, parametric resonance is 
rather ineffective for derivatively coupled scalar fields, no 
matter whether heavy or light. 



C. Preheating into Gauge Fields 

Our analysis so far indicates that in natural models 
of inflation the population of the universe with particles 
faces several challenges. First, if matter fields are heavier 
than the inflaton, parametric resonance is absent, and the 
inflaton cannot decay perturbatively into matter for kine- 
matical reasons. And second, if matter fields are lighter 
than the inflaton, parametric resonance or perturbative 
decay are possible, but then long-wavelength fluctuations 
of matter are produced during inflation. These conclu- 
sions crucially hinge on the assumption that matter fields 
can be represented by minimally coupled scalar fields, 
which, of course, is not the case. 

There are just two other types of matter fields the infla- 
ton can couple to: fermions (spin 1/2 and 3/2) and vec- 
tors (spin 1). Reheating into fermions is not very effec- 
tive because of the exclusion principle ^25] , so we are left 



1 



1 



1 



(32) 

Here = d^A^ — d^jA^ is the field strength of the 
gauge fields, A^i, = e'^''^'^ Ape /2 is its dual tensor, and F 
is a symmetry breaking scale, with dimensions of mass. 
Despite the appearances, the field (p couples derivatively 
to the gauge field, because A^^A^^'^ is a total derivative. 
Effective theories of this type appear in axion-like models 
P7] , in which the coupling of the scalar to the electromag- 
netic field arises from the triangle anomaly. This sort of 
coupling between a pseudoscalar and electromagnetism 
has been suggested as a mechanism to generate primor- 
dial magnetic fields in the early universe [2S], although 
a detailed study of coherent magnetic field production 
during inflation and the subsequent reheating stage has 
concluded that this generation is not sufficiently strong 
on the required scales [25] . 

To analyze the dynamics of reheating in this model, we 
choose to work in the transverse gauge (where diA^ — Q 
and A° = 0). In this gauge, the equations of motion for 
the inflaton and the non- vanishing components of are 



ml^+—e^='^MdjAu (33) 
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We expand Ai in terms of mode functions: 



(35) 



where cj^: = |k| = fc is the dispersion relation of the 
photons, and er(k) are circular polarization vectors, 
k X e± = =p i fc e± . Substituting this expansion into ( 34 ) 
yields the decoupled mode equations 



A, 



HAa 



k 4) 

aF 



A+ = 0. 



(36) 



As was the case in the previously considered mod- 
els, equation (36 1 exhibits parametric resonance for cer- 
tain values of k and F~^, which can be interpreted as 
explosive particle production (e.g., photon production). 
This can be seen analytically in the regime where back- 
reaction is neglected by putting ( [36| into the form of the 
Mathieu equation ([9|. The rescaled field variable, 

i± = A± (37) 

satisfies the Mathieu equation ^ with parameters 



(5 = 



/ 2k 



2 



and e — ± 



k 1 



mA,a t F ' 



(38) 



where we have introduced the dimensionless time vari- 
able t 2m^^t. There is a significant difference that dis- 
tinguishes couplings to gauge fields from those to scalar 
fields. Whereas for the latter the value of e is constant 
in a non-expanding universe, for the former it is propor- 
tional to the energy of the photons. Hence, e can become 
very large. Since the effective theory we are using breaks 
down at physical momentum k/a « F we find that at 
the end of infiation e can be as large as (t>c/m^ w 10^. 
However, this does not necessarily imply that paramet- 
ric resonance is very effective. In Fig. [8] we show the 
trajectories of different modes in the 5-e plane for dif- 
ferent values of k and (j)c/F — 1. As seen in the figure, 
although a mode can cross several instability bands, it 
does so in the narrow band regime, close to the tips of 
the instability bands, mainly because e rapidly decays as 
time evolves. 

In order to estimate the efficiency of preheating in 
this regime, let us estimate how long a mode remains 
in the first instability band (the one whose tip touches 
e = 0, S = 1), and how much its mode function grows 
during that time. For small e, the characteristic expo- 
nent along the first band is 



(39) 



Hence, its boundary is \S — 1\ = e, along which /x = 0. Let 
us write the time at which a mode enters and leaves the 



band as =F At, where ti is the time at which the mode 
crosses the center of the band, at 6 = 1. Using equations 
(38 1, and neglecting the factor 2/9t^, we find 



2k 



3/2 



At 



F 



(40) 



where ~ 1 is the time at the end of inflation, and Oe 
denotes the value of the scale factor at that time. 




FIG. 8: Instability diagram for solutions of the Mathieu equa- 
tion in the 5-e plane. The trajectories follow the evolution of 
the parameters in equation ( 38 1 starting at the end of infla- 
tion. The corresponding values of fee = k/{m^ac) are fee — 3/4 
(continuous), fee ~ 3/2 (dashed) and fee = 3 (dot-dashed). In 
all cases (pc/F — 1. 

We shall estimate now what is the growth of the mode 
function during that time. Inside the band, we can ne- 



glect the term {S— 1)^ in equation (39 1, so that fi ~ ei/2 



where ei is the value of e at time ti. Hence, we find that 
A^ grows by a factor 



exp(ei • At) « exp 



16ip F2 



2k 



3/2' 



(41) 



while the growth in A-j- is suppressed by an additional 
factor {tc/tY/'^. Note that the above expression is strictly 
valid only for At <IC ii, and provided the change in 
/i — e/2 at time ti is adiabatic. 



20 F 



< I. 



(42) 



Equation ( 41 1 hence implies that preheating into gauge 



fields is very effective for (t)c/F ^ 1. However, this is not 
the regime that applies in natural models of inflation, 
where (pc/ F is of order one (a small number for our pur- 
poses). In order to estimate the growth in the opposite 
regime, we shall set e oc (pc/F 0. In that case, the 
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FIG . 9: Num erical solution of the system of equations ([5|, 
(33 1 and (361 (without back-reaction) for k = m^. There is 
no amplification of the mode function, which is well-described 
by the approximate solution (|43[). 



solution of the Mathieu equation has the growing mode 



6kt 



2/3 



(43) 



which implies that the original variable oscillates with 
constant amplitude. Therefore, there is no parametric 
resonance in the decay of the inflaton into gauge fields. 
As before, the lack of parametric resonance is due to the 
modest value of (l>c/F, which is a reflection of the weak 
couplings between the inflaton and matter. Numerical 
investigation including the effects of expansion confirms 
this. Fig. ([9| shows for instance the evolution of yl+ for 
k ~ which is very well approximated by equation 
(43 1. As clearly seen in the figure, there is no particle 



production. 

Because the equation of motion is invariant under t — > 
t + TT and e — > — e, similar results hold for A_ modes. 
Notice that the comoving number density is now given 

by 



a 

2uji, 



(44) 



IV. CONCLUSIONS 

It is obviously necessary to repopulate the universe 
with matter after a period of early universe inflation. 
The detailed process through which this takes place can 
be of great significance since, although a thermally equi- 
librated bath ultimately results, non-equilibrium relics 
from preheating can have important cosmological conse- 
quences, even if the ultimate reheat temperature of the 
universe is below the mass scale of any relics of interest. 



The possible implications of preheating include the 
overclosure of the universe through the production of 
monopoles, moduli or gravitinos, and new possibilities 
for the generation of the baryon asymmetry of the uni- 
verse, either at the Grand Unified (GUT) [30^,31^ or elec- 
troweak [32l [33l [Ml |35] scale. It is therefore important 
to examine the different ways in which parametric reso- 
nance may take place in inflationary models. 

Because phenomenologically viable inflation requires 
extremely flat potentials, models in which the couplings 
of the inflaton to itself and to other flelds are natu- 
rally suppressed are particularly attractive. One way 
to achieve this is to construct the inflaton as a pseudo- 
Nambu-Goldstone boson of a spontaneously broken U{1) 
symmetry — natural inflation. The U{1) symmetry is 
then realized as a shift symmetry on the inflaton field, 
and explicit soft breaking terms render this symmetry 
approximate, and generate a naturally small mass and 
approximately flat potential for the inflaton. A conse- 
quence of the residual approximate symmetry is that di- 
rect couplings o the inflaton to matter flelds are corre- 
spondingly suppressed, leading to the interesting possi- 
bility that derivative couplings, unconstrained by a shift 
symmetry, may be the dominant type of coupling to mat- 
ter. 

In this paper, we have used analytic and numerical 
techniques to study reheating in models in which deriva- 
tive couplings between the inflaton and matter fields are 
expected to play the dominant role — such as in natu- 
ral infiation. We have seen that successful reheating 
places non-trivial constraints on these models and con- 
fronts them with several challenges. In particular, for 
heavy scalar fields (to^ ^ m^) parametric resonance is 
ineffective in producing matter particles, and the pertur- 
bative decay is forbidden on kinematic grounds. For light 
scalars {rrix ^ there is a a novel long- wavelength in- 
stability that causes the zero mode of matter to grow as 
a power law, but in this case, the production of fluctu- 
ations of matter during inflation leads to several cosmo- 
logical problems. We note that we have not addressed 
any mechanism to keep matter fields light, which could 
lead to further constraints. 

Some of these problems are avoided when the infiaton 
couples derivatively to gauge fields, as in axion-like mod- 
els. Because of the gauge symmetry, these vector fields 
are automatically light, thus allowing the perturbative 
decay of the inflaton. At the same time, the confor- 
mal nature of their couplings to gravity circumvents all 
the problems associated with the presence of light scalars 
during inflation. Nevertheless, we have found that para- 
metric resonance is absent for the values of parameters 
implied by natural inflationary models. As for couplings 
to scalars, the origin of this absence is the weakness of 
the coupling of the inflaton to matter, which is charac- 
terized by the ratio of the inflaton at the end of inflation 
to the spontaneous breaking scale, and is hence of order 
one at most. 
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